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Abstract. We report the results of an XMM-Newton observation of the Lagoon Nebula (M 8). Our EPIC images of thrs region 
reveal a cluster of point sources, most of which have optical counterparts inside the very young open cluster NGC 6530. The 
bulk of these X-ray sources are probably associated with low and intermediate mass pre-mam sequence stars. One of the sources 
experiences a flare-like increase of its X-ray flux making it the second brightest source in M 8 after the 04 star 9 Sgr. The X-ray 
spectra of most of the brightest sources can be fitted with thermal plasma models with temperatures of kT ~ a few keV. Only 
a few of the X-ray selected PMS candidates are known to display Ha emission and were previously classified as classical 
TTauri stars. This suggests that most of the X-ray emitting PMS stars in NGC 6530 are weak hne T Taun stars, n a ltion to 
9 Sgr our EPIC field of view contains also a few early-type stars. The X-ray emission from HD 164816 is found to be typical 
for an 09.5III-IV star. At least one of the known Herbig Be stars in NGC 6530 (LkHa 115) exhibits a relatively strong X-ray 
emission, while most of the main sequence stars of spectral type B 1 and later are not detected. We also detect (probably) diffuse 
X-ray emission from the Hourglass Region that might reveal a hot bubble blown by the stellar wind of Herschel 36, the ionizing 
star of the Hourglass Region. 
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1. Introduction 

The Lagoon Nebula (M 8 = NGC 6523) is one of the closest 
and brightest H n regions in the Galaxy. The nebula is ionized 
by the 04 V star 9 Sgr and the binary system HD 165052 
(06.5 V + 06.5 V). M 8 harbors also the very young open clus- 
ter NGC 6530 as well as a compact H 1 1 region, the so-called 
Hourglass Region (hereafter HG), which is ionized by the 07 V 
star Herschel 36. NGC 6530 is believed to be at the origin of 
a sequential star formation process inside the Lagoon Nebula 
(Lightfoot et al. 1984). NGC 6530 most probably formed from 
a pre-existing massive molecular cloud. The ionizing radiation 
of the newly formed hot stars in the cluster then created a 
cavity allowing us to see inside the parent molecular cloud. 
Furthermore, NGC 6530 may have triggered the formation of 
9 Sgr, which subsequently caused the formation of the massive 
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protostar Herschel 36 (= HD 164740, hereafter H 36). 

According to van den Ancker et al. (1997), star formation 
within NGC 6530 started a few 10 7 years ago and is probably 
still ongoing. Sung et al. (2000) found that 9 Sgr and H 36 are 
well fitted to an isochrone of age ~ 1.5 Myr, while the age of 
the stars at the fainter end of the main-sequence of NGC 6530 
is about 4 Myr and the massive PMS star LkHa 112 should be 
only about 0.4 Myr old. Using UBVRI and Ha photometry of 
NGC 6530, Sung et al. (2000) identified 58 pre-main sequence 
(PMS) stars with Ha emission as well as 17 PMS candidates. 

Many recent studies of the Lagoon Nebula focused on the 
Hourglass Region as seen in different wavebands. The HG 
nebula is believed to be an ionized cavity in an inhomogeneous 
clumpy molecular cloud. The apparent bipolar structure in the 
optical band most likely results from a non-uniform extinction 
along the line of sight. Woodward et al. (1986) proposed a 
model in which two high density knots on the eastern side 
protect the material behind them from the ionizing radiation 
while the western edge of the HG would be defined by a 
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foreground cloud obscuring H 36. 

In the present paper, we discuss an XMM-Newton observa- 
tion of M 8, that reveals a number of point sources towards the 
open cluster NGC 6530. In Sects. 2 and 3.1, we investigate the 
light curves and spectra of the brightest sources and we dis- 
cuss their possible association with PMS stars in NGC 6530. In 
Section 4, we discuss also the soft X-ray emission from the HG 
nebula. Finally, our conclusions are highlighted in Section 5. 

2. X-ray sources inside the Lagoon Nebula 

The Lagoon Nebula was observed serendipitously for about 
20ksec with the XMM-Newton observatory (Jansen et 
al. 2001) during our AOl observation of the 04 V star 9 Sgr 
(JD 245 1976.978- 2451977.239). The two EPIC-MOS instru- 
ments were operated in the full frame mode (Turner et al. 2001) 
whilst the EPIC-pn camera was used in the extended full frame 
mode (Struder et al.2001). All three EPIC instruments used 
the thick filter to reject optical light. We used version 5. 1 of the 
XMM-Newton Science Analysis System (SAS) to reduce the 
raw EPIC data. More details on the pipeline processing of the 
data are given in the paper dealing with the 9 Sgr data (Rauw 
et al. 2002). 


The combined EPIC soft band (0.5 - 1.2 keV) X-ray 
image around 9 Sgr reveals a number of point-like (and maybe 
also diffuse) sources (see Fig. 1). Figure 2 displays the EPIC 
contours of the X-ray sources superimposed on a Digital Sky 
Survey image of M 8. This image indicates that most of the 
fainter sources are located in the region of the very young 
open cluster NGC 6530. One rather bright and apparently 
diffuse source (to the south-west of 9 Sgr) is associated with 
the Hourglass Region. 

Because of the contamination of the EPIC data by hard 
stray light photons from the nearby LMXB GX 9 + 1 , the source 
detection was performed on the soft band images (i.e. where 
there is no contamination) only. For details on the properties of 
the straylight, we refer again to the work of Rauw et al. (2002). 
Using the SAS source detection algorithms (see e.g. Hasinger et 
al. 2001), we find 250 sources that are detected simultaneously 
in the MOS 1 , MOS2 and pn images with a combined likelihood 

3 

~ > 9.0 

t“l 

where the summation runs over the MOS1, MOS2 and pn im- 
ages. In this formula, pi stands for the probability that a random 
Poissonian fluctuation may have caused the observed source 
counts in the image of the detector i. We then inspected each 
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Fig. 2. EPIC X-ray contours (corresponding to the combined MOS1 + MOS2 + pn data in the 0.5 1.2 keV range) superimposed on a DSS 

optical image of the Lagoon Nebula. The different contours correspond to 5, 15, 24, 50 and 200 counts in the X-ray image. 


source manually to reject spurious detections due e.g. to hot 
pixels. This resulted in a total of 220 confirmed detections. Out 
of these 220 sources, 1 17 have a combined likelihood > 20.0 
and are detected in the individual images from all three EPIC 
instruments with — lnpj > 3.0. To this list we add two sources 
that are not detected in at least one of the instruments because 
they fall outside its field of view, but are clearly detected in 
the remaining instruments. This gives us a total of 1 19 highly 
significant detections. In the following, we will restrict our 
discussion to these latter objects. The faintest sources in this 
category have about 10” 3 cts s -1 over the 0.5 - 1 .2 keV band 
of the MOS instruments. Assuming a 1 keV thermal spectrum 
with a neutral hydrogen column density of 0.17 10 22 cm 
(see also Sect. 3.1 below), the faintest sources correspond to 
an observed flux of about 8.9 10~ 15 ergcm 2 s 1 and an un- 
absorbed flux of 13.1 10~ 15 erg cm -2 s" 1 in the 0.5 - 5.0keV 
energy range. Note that the corresponding observed flux in the 
0.5 - 1.2 keV soft band would be 5.2 10 -15 erg cm -2 s -1 . 

We have cross-correlated the positions of the 1 19 sources 
with the results of the UBVR1 and Ha photometry of NGC 6530 
presented by Sung et al. (2000). We searched the Sung et al. 
catalogue and the SIMBAD database for stars within a radius 
of 12" (corresponding roughly to twice the FWHM of the on- 
axis point spread function (PSF) of the XMM mirror modules 
at 1 .5 keV; Jansen et al. 2001) from the coordinates of the X-ray 
sources: 

• the brightest X-ray source is associated with the 04 V((f + )) 
star 9Sgr. The data of this star have been analysed else- 
where (Rauw et al. 2002) and we will not repeat this dis- 
cussion here. 

• a probably diffuse X-ray source is associated with the 
Hourglass Region. We shall return to the HG in Sect. 4. 


• 63 X-ray sources have a single optical counterpart within 
a radius of less than 12". The average angular separation 
between the X-ray source and the optical counterpart is 
(4.1 ± 2.4)". The brightest soft X-ray sources in this cate- 
gory are HD 164816 (09.5 III-IV) and SCB 731. The latter 
is a very red 16th magnitude star without outstanding fea- 
tures (such as Ha emission) in its optical photometry. The 
X-ray source associated with SCB 731 appears elongated 
since it lies rather far off-axis. 

• 30 X-ray sources have two or more optical stars falling 
within the 12" radius. Thirteen of these sources have at 
least one known Ha emission star inside their ‘error box’. 

• 24 sources do not seem to have a known optical counterpart 
in the NGC 6530 cluster. 

The sources of the last three categories are listed in 
Tables 1, 2 and 3 respectively. We have adopted the naming 
conventions recommended by the XMM SOC and the IAU (see 
the XMM-Newton Newsletter #4, April 2001): the XMMUJ 
prefix is followed by the right ascension HHMMSS.s (in hours, 
minutes, seconds and tenths of seconds) and the declination 
of the source +/-DDMMSS (in degrees, arcminutes and 
arcseconds). Note that the coordinates refer to equinox 2000.0 
and are truncated, not rounded. 

We have extracted light curves and spectra for the bright- 
est objects that are not affected by the straylight. Except for 
sources near the edge of the field of view, we used a circular 
extraction region with a radius between 19 and 40" depending 
on the angular separation between neighbouring X-ray sources. 
For sources near the edge of the field of view, we used ellipti- 
cal extraction regions to match as closely as possible the shape 
of the point spread function. The light curves and spectra were 
accumulated over the PI range 200 to 10000 (corresponding 
roughly to E e [0.2, 10.0] keV), except for objects that are 
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Table 1. X-ray sources detected in NGC 6530 and found to have a single counterpart in the catalogue of Sung et al. (2000) or the SIMBAD 
database The first column yields the name of the X-ray sources according to the naming conventions for serendipitous XMM-Newton sources 
columns [2] and [3] provide the vignetting and background corrected EPIC-MOS countrates in the 0.5 - 1.2 keV band. The number of the optical 
counterpart from Sung et al. (2000) and the separation from the position of the X-ray source are given in columns [4] and [5] The photometric 
data are taken from Sung et al. (2000), except for V 1752 Sgr. Column [8] indicates whether variability ,s detected at the 95% confidence level 
(Y) or not (N). Those sources flagged with a T show some hints of variability in their light curve, but this is not found to be significant bv 
the Kolmogorov-Smimov test. Finally, the last column yields information about the spectral types, pre-main sequence classification and the 
membership probability (P^) according to the proper motion study of van Altena & Jones (1972). 


XMMUJ 

MOS1 

MOS2 

SCB 

d 

V 

B - V 

Variability 

Comments 

[i] 

10 3 ctss -1 

10~ 3 ctss -1 


(arcmin) 





[2] 

[3] 

[41 

[5] 

[6] 

[7] 

[8] 

[9] 

180253.7-242053 

4.1 ± 1.0 

2.8 ± 1.1 


0.14 

15.80 



V 1752 Sgr= PMS 

180322.8-242208 

3.0 ± 0.6 


7 

0.07 

10.62 

0.96 

N 

180323.3-241525 

1.4 ±0.7 

1.5 ±0.7 

20 

0.12 

15.17 

1.28 


180326.9-241246 

2.9 ± 0.7 

2.1 ±0.6 

47 

0.02 

11.83 

0.26 


= 0.80 

180337.5-242116 

0.5 ± 0.4 

1.1 ±0.4 

1001 

0.17 

17.32 

1.31 



180340.1-241923 

1.6 ±0.4 

1.7 ±0.4 

138 

0.05 

15.05 

1.24 



180341.5-241222 

1.5 ±0.5 

1.0 ±0.4 

151 

0.03 

15.38 

1.33 



180341.5-242416 

1.9 ±0.4 

2.2 ± 0.4 

157 

0.09 

15.44 

1.47 


PMS 

180343.4-242340 

2.6 ±0.5 

2.0 ± 0.4 

174 

0.06 

16.90 

1.44 


180345.1-242205 

5.1 ±0.6 

7.3 ± 0.7 

182 

0.00 

9.95 

0.18 

7 

B2.5V, Pp. = 0.86 

180347.0-241647 


1.6 ±0.4 

199 

0.04 

10.75 

0.09 



180347.2-242201 

1.0 ±0.6 

1.8 ± 0.7 

202 

0.06 

16.64 

1.55 



180349.1-241537 

0.7 ± 0.3 

0.8 ± 0.3 

220 

0.06 

16.87 

1.31 



180349.2-242343 

0.8 ± 0.4 

0.7 ±0.3 

223 

0.11 

16.85 

1.35 



180350.9-242018 

2.5 ±0.5 

1.7 ±0.4 

238 

0.10 

16.96 

1.44 



180352.6-242256 

4.3 ± 0.6 

4.6 ± 0.6 

253 

0.09 

13.80 

1.23 

N 

P„ - 0.00 

180355.4-241701 

2.8 ± 0.5 

4.7 ± 0.6 

274 

0.04 

14.61 

1.06 

N 


180356.6-241845 

14.3 ± 1.0 

16.6 ± 1.1 

282 

0.06 

6.99 * 

0.02 

N 

HD 164816 (09.5 II1-IV) 

180357.3-241611 

1.3 ±0.4 

1.3 ±0.4 

285 

0.03 

16.27 

1.20 


180358.3-242455 

3.7 ±0.5 

3.4 ± 0.5 

295 

0.09 

15.51 

1.40 

N 


180358.6-242530 

6.4 ± 0.7 

6.5 ± 0.7 

290 

0.17 

12.64 

1.03 

7 

SV Sgr = PMS, P^ =0 00 

180400.2-241915 

3.2 ± 0.5 

4.0 ± 0.6 

314 

0.07 

15.28 

1.14 

N ■ 


180402.9-242140 

4.4 ± 0.6 

5.4 ± 0.6 

340 

0.01 

14.85 

1.35 

N 


180404.9-241643 

1.3 ±0.4 

0.7 ± 0.3 

360 

0.07 

16.31 

1.40 


180407.2-241926 


1.3 ±0.4 

379 

0.05 

10.76 

0.25 


P» = 0.00 

180407.9-242126 

1.7 ±0.4 

2.0 ± 0.4 

387 

0.00 

16.43 

1.32 



180409.3-242017 

1.7 ±0.5 

1.4 ±0.4 

412 

0.19 

11.60 

0.25 


P» = 0.62 

180409.4-242715 

2.4 ±0.5 


410 

0.12 

16.69 

1.60 



180411.0-242149 

2.0 ±0.4 

1.6 ±0.4 

426 

0.08 

10.37 

0.12 


B3Ve, P^ = 0.84 

180411.8-242527 

1.7 ±0.4 

2.6 ± 0.5 

432 

0.06 

11.69 

0.24 


A0, P^ = 0.52 

180413.0-241825 

1.8 ±0.4 

1.7 ±0.4 

450 

0.06 

14.34 

0.98 



180413.2-242452 

1.8 ±0.4 

1.8 ±0.4 

452 

0.05 

15.31 

1.35 



180413.8-242619 

1.4 ±0.5 

1.3 ±0.5 

457 

0.08 

16.32 

1.53 



180413.9-242405 

2.6 ± 0.5 

2.6 ±0.5 

454 

0.03 

15.59 

1.40 



180414.3-241938 

0.8 ± 0.5 

1.3 ±0.5 

465 

0.05 

14.39 

1.22 



180414.4-242155 

7.6 ± 0.7 

7.2 ± 0.8 

470 

0.07 

13.45 

1.20 

N 

Py. - 0.00 

180415.1-242112 

1.6 ±0.4 

2.5 ± 0.5 

480 

0.13 

16.83 

1.49 



180415.8-241457 

1.9 ±0.5 

1.3 ±0.5 

481 

0.04 

12.79 

0.27 


Pn = 0.76 

180416.8-242835 

180417.5-242415 

2.4 ± 0.6 
2.0 ± 0.5 

1.7 ±0.5 
3.2 ± 0.5 

502 

504 

0.03 

0.12 

11.56 

14.26 

0.23 

1.22 


A0, = 0.60 

180420.4-242204 

1.2 ±0.5 

1.6 ±0.6 

546 

0.06 

15.43 

1.33 



180420.9-242324 

2.7 ± 0.5 

4.3 ± 0.6 

556 

0.07 

16.98 

1.41 

N 

PMS“ 

180421.1-241951 

4.2 ± 0.6 

4.1 ±0.6 

557 

0.04 

14.71 

1.11 

N 

180422.6-242355 

5.2 ± 0.7 

6.2 ±0.8 

586 

0.08 

16.81 

1.47 

N 


180423.6-242951 

1.4 ±0.4 

0.8 ± 0.5 

592 

0.06 

16.76 

1.42 



180424.2-241829 

3.2 ± 0.6 

1.7 ±0.5 

600 

0.03 

14.19 

1.04 



180424.5-242025 

2.6 ± 0.5 

2.4 ±0.5 

607 

0.05 

16.47 

1.45 



180425.8-241750 

1.3 ±0.5 

1.3 ±0.5 

617 

0.05 

15.36 

1.22 




Another star. Walker 54 (V = 11.6) lies at 0.18' according to van den Ancker et al. (1997). 
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XMMU J 
[1] 

MOS1 

10~ 3 cts s -1 
[2] 

MOS2 
10 -3 ctss -1 

[3] 

SCB 

[4] 

d 

(arcmin) 

[5] 

V 

[6] 

B-V 

[7] 

Variability 

[8] 

Comments 

[9] 

180425.8-242357 

2.8 ± 0.8 

2.3 ± 0.7 

627 

0.14 

16.69 

1.43 



180428.6-241855 


3.9 ± 0.7 

654 

0.04 

14.17 

0.97 

N 


180428.7-241400 

1.1 ±0.6 

2.6 ± 0.8 

656 

0.06 

11.99 

0.52 


= 0.83 

180428.7-241444 

1.7 ±0.5 

1.2 ±0.5 

655 

0.05 

15.39 

0.63 



180429.5-242021 

3.7 ± 0.6 

2.8 ± 0.6 

665 

0.11 

16.86 

1.54 

N 


180431.9-242406 

2.6 ± 0.6 

2.7 ± 0.6 

700 

0.11 

14.59 

1.18 


5 

180432.4-242122 

2.9 ± 0.6 

2.5 ± 0.6 

698 

0.11 

16.11 

1.41 



180432.9-241844 

3.4 ± 0.7 

3.1 ±0.6 

707 

0.03 

10.42 

* 0.12 

N 

B2.5Vne, P ^ = 0.75 

180434.3-241925 

2.2 ± 0.6 

3.0 ±0.6 

725 

0.04 

12.85 

0.87 


P ^ = 0.01 

180435.5-242729 

15.5 ± 1.5 

15.7 ± 1.5 

731 

0.03 

16.53 

1.50 

Y 


180437.1-242956 

2.2 ± 0.6 

1.2 ±0.6 

747 

0.03 

15.02 

1.40 



180438.1-242043 

1.9 ±0.6 

2.1 ±0.6 

761 

0.07 

15.41 

1.40 



180439.1-242409 

1.4 ±0.7 

2.3 ± 0.7 

763 

0.08 

10.21 

1.21 


K0III, Pp = 0 00 

180441.8-241855 


2.6 ± 0.7 

786 

0.06 

14.62 

1.04 



180450.6-242540 

4.3 ± 1.0 


879 

0.03 

11.94 

0.39 


LkHa 115 (B4Ve), P^ = 0.79 


6 Another star. Walker 82 ( V = 12.4) lies at 0. 1 8' according to van den Ancker et al. (1997). 

Table 2. Same as Table 1 for X-ray sources with two or more counterparts from Sung et al. (2000) inside a radius of 1 2 around the position of 
the X-ray source. 


XMMUJ 

[1] 

MOS1 

10 -3 ctss -1 

[2] 

MOS2 
10 -3 ctss -1 

[3] 

SCB 

[4] 

d 

(arcmin) 

[5] 

Variability 

[6] 

Comments 

[7] 

180323.8-241556 

2.7 ± 0.6 

3.1 ±0.6 

28/14 

0.10/0.16 



180338.1-242235 

1.8 ±0.5 

1.3 ±0.5 

1039/1028/1036/1027 

0.03/0.06/0.18/0.19 



180339.8-242200 

1.6 ±0.4 

1.7 ±0.4 

1032/1033 

0.15/0.19 


1032, 1033 = PMS 

180340.3-242316 

4.9 ± 0.6 


142/146 

0.07/0.09 


142 ~ LkHa 106 

180350.5-242110 



240/228/229 



240 = LkHal08 b 

180407.4-242226 

7.8 ± 0.8 

7.8 ± 0.8 

382/380 

0.05/0.07 

N 


180407.7-24231 1 

1.7 ± 0.4 

1.4 ±0.4 

388/390 

0.05/0.05 


390 = PMS 

180412.5-241950 

4.0 ± 0.6 

2.9 ± 0.5 

440/436 

0.13/0.16 

N 

440 = PMS 

180414.4-242332 

2.2 ± 0.6 

1.8 ±0.6 

458/476 

0.12/0.18 



180414.7-242224 

1.5 ±0.6 

3.2 ± 0.6 

468/1005 

0.07/0.07 


1005 ~ PMS 

180415.3-241902 

0.8 ± 0.4 

1.7 ±0.5 

478/482 

0.08/0.11 


482 = PMSc 

180420.2-242250 

2.7 ± 0.5 

2.8 ± 0.5 

540/531 

0.04/0.19 


540 = PMS, 531 = PMSc 

180420.4-242810 

2.3 ± 0.6 

2.4 ± 0.6 

552/547 

0.15/0.16 


552, 547 = PMS 

180421.8-242114 

1.4 ±0.8 

1.8 ±0.8 

573/561/558 

0.09/0.13/0.14 



180421.8-242039 

4.2 ± 0.6 

4.7 ± 0.7 

571/572 

0.11/0.18 

N 


180423.0-242116 

7.5 ± 0.8 

8.3 ± 0.9 

582/593/587 

0.12/0.13/0.14 

N 


180425.0-242747 

2.6 ± 0.6 


611/618 

0.04/0.18 



180425.7-242254 

1.9 ± 0.7 

3.7 ± 0.9 

614/625 

0.11/0.19 


625 = LkHa 1 13 

180426.9-242100 

2.0 ± 0.5 

2.3 ±0.5 

631/635 

0.11/0.11 



180427.5-242354 

1.9 ±0.6 

2.0 ± 0.6 

641/632/648 

0.03/0.17/0.18 



180429.0-242225 

5.1 ± 0.7 

6.5 ± 0.8 

658/672/673 

0.12/0.15/0.16 

N 


180429.1-242147 

1.5 ±0.6 

2.0 ± 0.7 

663/657 

0.09/0.10 



180429.5-242300 

7.7 ± 0.9 

7.0 ± 0.9 

671/659 

0.06/0.13 

N 


180430.7-242634 

4.6 ± 0.8 

3.7 ± 0.7 

681/682 

0.04/0.06 

N 


180433.6-242706 

6.8 ± 1.0 

3.5 ± 0.7 

718/713 

0.07/0.21 

N 

713 = LkHa 114 

180434.2-242154 

1.7 ± 0.5 

1.6 ±0.6 

723/720/716 

0.08/0.11/0.14 



180436.0-241959 

1.0 ±0.6 

2.3 ± 0.7 

740/738 

0.05/0.13 



180439.0-241924 

2.0 ± 0.6 

2.0 ± 0.6 

765/773 

0.03/0.14 



180441.8-242050 

1.4 ±0.5 


787/790 

0.10/0.15 



180448.7-242634 

10.5 ± 1.4 

9.3 ± 1.4 

864/862/871 

0.02/0.11/0.19 

N 

862 = PMS 


a SCB 1027 & 1028 are classified as PMS stars, while SCB 1036 & 1039 are PMS candidates (Sung et al. 2000). 

* XMMU J1 80350.5— 2421 10 lies within the wings of the PSF of 9 Sgr. Therefore no accurate count rate can be provided. 
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contaminated by the stray light. In the latter cases, we extracted Table 3. X-ray sources detected in NGC 6530 that have no counterpart 

light curves only (i.e. no spectra) over the PI range 200 to 1 500 in the Sung et al. (2000) catalogue and the SIMBAD database within 

(i.e. E € [0.2, 1.5] keV). a radius of 12". 


2. 1. Light curves 

For all the light curves we adopted time bins of 400 s. The 
light curves were tested for variability using a Kolmogorov- 
Smirnov test against a model of constant count rate. Among 
the sources with a single optical counterpart, SCB731, the 
second brightest X-ray source, is clearly variable: the light 
curve (extracted over an elliptical region adapted to the shape 
of the PSF) reveals a strong X-ray flare towards the end of 
our observation (Fig. 3). The X-ray flux suddenly rises by 
nearly two orders of magnitude and after about 5000 s starts to 
decline again. 

SCB 182 and SCB290 show some marginal variability in 
their light curves. The behaviour is the same in all three EPIC 
instruments, but the variability occurs at a rather low level 
so that it is found not to be significant by the Kolmogorov- 
Smirnov test (at the P > 95% confidence level). 

None of the other relatively bright sources turned out to display 
significant variability (95% level) over the duration of our ob- 
servation. 

We reduced also a ROSAT - PSPC observation of M8 
(rp900374n00, integration time 10.4 ksec, JD 2449078.532 - 
2449079.569) using the XSELECT software. In this exposure, 
SCB 253 appears as the third brightest source (after 9 Sgr and 
HD 164816) with a net count rate of (6.9 ± 1.0) 10~ 3 ctss _1 
over the entire energy range of the PSPC instrument. Therefore 
it seems likely that SCB 253 experienced a flare at the time 
of the ROSAT observation. Other sources that are detected 
in this PSPC image are the Hourglass Region, SCB 7, 47, 
182, 274, 340, 452/454, 470, XMMU J180323.8-241556, 

180407.4-242226, 180253.3-241955, 180356.4-242340,’ 

180339.1—243302 as well as a source which is not seen with 
XMM at RA = 18:04:03.2, DEC = -24:31:08.5 (equinox 
2000 . 0 ). 

2.2. Spectra 

For the spectra, we adopted the redistribution matrices pro- 
vided by the EPIC instrument teams (versions available in May 
2001) and we used the SAS to build the appropriate ancillary 
response files for each EPIC instrument. The spectra were 
binned to reach a minimum of 10 counts per channel and the 
background corrected spectra were analyzed using the XSPEC 
software (version 11.00). Because of the strong noise in the 
pn detector below 0.2 keV and the uncertain calibration of 
the EPIC instruments at low energies, we ignored the binned 
energy channels below 0.2 ke V for the XSPEC fits. Note that for 
the X-ray sources associated with HD 164816 and SCB 731, 
we have a total of respectively 300 and 800 net counts for each 
MOS spectrum. The pn spectra contain about twice as much 


XMMU J MOS1 MOS2 Variability 

10~ 3 ctss _1 10~ 3 ctss _1 

W [2] [3] [4] 

180252.9-241525 3 ~ ' 

180253.3 -241955 6.4 ±1.4 N 

180257.3- 242014 3.9 ±1.1 2.9 ±1.1 

180339.1- 243302 5.1 ± 10 6.9 ±1.1 N 

180351.3- 241149 1.8 ±0.5 1.8 ±0.6 

180356.4- 242340 1.3 ± 0.4 2.3 ± 0.4 

180356.5- 242018 0.8 ±0.4 0.4 ± 0.4 

180357.1- 242055 1.8 ±0.4 

180402.6- 242645 2.3 ± 0.5 2.2 ± 0.4 

180404.5- 241941 0.7 ±0.3 1.2 ± 0.4 

180405.0- 242511 2.9 ±0.5 3.4 ± 0.5 

1 80405.2- 242340 0.8 ± 0.3 1 . 1 ± 0.4 

180409.3- 241622 1.2 ±0.4 0.8 ± 0.4 

180412.2- 241907 1.5 ±0.5 1.3 ±0.5 

180413.2- 241517 1.4 ±0.5 1.0 ±0.5 

180413.9- 241722 0.8 ± 0.4 1.3 ±0.5 

1 804 1 9.6- 242603 1 . 1 ± 0.4 1 . 1 ± 0.4 

180420.3- 241920 0.9 ± 0.4 1.6 ±0.5 

180424.7- 242454 1.9 ± 0.6 1.5 ± 0.5 

180430.5-242850 1.2 ±0.5 1.9 ±0.6 

180435.0- 242334 3.3 ± 0.7 4.2 ± 0.7 N 

180440.3- 242643 7.3 ±1.1 6.9 ±1.1 N 

180441.9- 242248 2.0 ±0.7 1.9 ±0.7 

180441.9- 242128 3.8 ± 0.9 1,7 ± 0,7 

° XMMU J 180252.9— 241525 lies outside the field of view of the 
MOS instruments, but is detected as a rather bright source with 
(1.14 ± 0.23) 10 -2 ctss -1 in the pn image. 


counts. 

We fitted the spectra using either an absorbed one- 
temperature mekal thermal plasma model (Mewe et al. 1985, 
Kaastra 1992) or an absorbed power law model. The best 
fitting model parameters are listed in Table 4 below. In several 
cases, the fits with a power law with a rather large photon index 
T are slightly better than those with the thermal plasma model. 
More complex models, such as two-temperature mekal or 
mekal + power law models could also improve the For 
instance, the spectrum of XMMU J 1 80407 .4-242226 clearly 
requires a 2-T thermal model. However, we caution that the 
quality of our spectra is usually not sufficient to distinguish a 
very hot (kT ~ a few keV) multi- temperature thermal plasma 
emission from a power law spectrum. 

For the 1-T mekal fits, we are usually left with a number of 
nearly equivalent solutions corresponding to either intrinsically 
harder, but less absorbed or more absorbed but intrinsically 
softer plasma models. Apart from the late O-star HD 164816 
(XMMU J1 80356.6— 241845), all the sources have rather hard 
spectra with kT > 0.65. 
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Table 4. Best fitting parameters for absorbed thermal plasma mekal and absorbed power law models fitted to the EPIC spectra of the brightest 
X-ray sources in NGC 6530. The observed fluxes quoted in columns [5] and [9] are evaluated over the energy range 0.5 - 5 keV. The quoted 
uncertainties correspond to formal errors on the spectral fits. 


wabs*mekal 


XMMUJ 

N h 

kT 

xi (d.o.f.) 

flux 


10 22 cm -2 

keV 


erg cm - 2 s 

[1] 

[2] 

[3] 

[4] 

[5] 

180345.1-242205 

0.22t;$ 

3.18t;f® 

0.98 (60) 

1.0910“ 

180356.6-241845 

0.50l°? 

0.21?;;$ 

1.09(152) 

0.8410“ 

180358.6-242530 

0.291;?? 

3.59l l 6 $ 

1.08 (48) 

2.22 10“ 

180414.4-242155 

0.l6t;°? 

3.05t;?? 

1.30(55) 

1.09 10“ 

180422.6-242355 

0.15±;° o ? 

2.95t;?? 

1.02(61) 

1.20 10“ 

180435.5-242729 

0.26±;$ 

6.64l 2 ;$ 

0.82(130) 

5.21 10“ 

180407.4-242226 

o.i7±;g 

/ 0.841;?? 1 
1 5.551?;?? J 

0.99 (63) 

1.08 10“ : 

180423.0-242116 

0.95±i| 

0.651;?? 

1.54 (66) 

0.89 10“ 

180429.0-242225 

i.06t;l! 

0.701; 2 ? 

1.55(40) 

0.6510“ 

180429.5-242300 

o.i5t;$ 

2.921;?? 

1.20 (90) 

1.1310" 

180448.7-242634 

i-04t;}? 

i-oii;?? 

1.51(45) 

1.2610“ 

180339.1-243302 

0-141;?? 

1.781;?? 

1.07 (36) 

0.6710“ 

180440.3-242643 

0.18t;$ 

3.491;?? 

0.91(106) 

1.71 10“ 


wabs* power 

T xi fd.o.f") flux 

erg cm - 2 s“ J 

[7] [8] [9] 


N h 

10 22 cm -2 

[6] 


0.43l$ 
0.43l;l 2 
0 341$ 

0.32l;$ 

0.4ll;J 2 

0.42l.j? 

0.35l;$ 

o.42!;13 

0.72l$ 

0.40-3 


2.781:31 

2.921;?? 

i.8ol:iS 

2.66l;33 

3.i5l;g 

2.991;?? 

2.62l;3| 

2.801:3 

4 . 8711:3 

2.65111 


1.07(55) 

0.82(61) 

0.78(130) 

1.35(65) 

0.95 (66) 
1.01(40) 
1.21(90) 
0.75 (45) 

1.24(36) 

0.79(106) 


1.08 10“ 13 

1.13 10“ 13 
5.27 10“ 13 

1.05 10' 13 

1.19 10“ 13 
0.83 10“ 13 

1.14 10“ 13 
1.58 10' 13 

0.55 10“ 13 
1.67 10“ 13 


The spectrum of HD 164816 could be better fitted by 
adding a second (harder) component to the model, but 
the properties of this second component are not well con- 
strained by the data. We note that the best fit model yields 
N h = 0.50 10 22 cm -2 which is significantly larger than the 
value of the interstellar column density 0.15;!; $ 10 22 cm 2 de- 
rived by Diplas & Savage (1994) from the interstellar Lya line. 
This suggests that part of the fitted absorption could actually 
arise in circumstellar (i.e. stellar wind) material. Correcting 
the X-ray flux for the pure interstellar column density only 1 , 
we infer an ‘intrinsic’ (i.e. at the top of the stellar wind) 
X-ray luminosity of 5.7 10 31 ergs -1 . Adopting a bolometric 
correction of -3.00, we infer logLx/Lbol = -6 96 which is 
in pretty good agreement with the empirical Lx - 1-bot relation 
for O-stars proposed by Berghofer et al. (1997). It is worth 
pointing out that Howarth et al. (1997) suggested that this star 
could be an SB2 binary. Provided that the interstellar column 
alone properly accounts for the low energy absorption, we find 
that the Lx/Lboi ratio yields no evidence for an excess X-ray 
emission attributable to a wind interaction in a colliding wind 
binary. 

If we assume that Nh = Nh,ism = 0.17 10 22 cm 2 (see 
Sect 3.1 below) for the other sources, we find that most fits 
yield kT of a few keV. Note that the spectrum of SCB731 


1 If we correct for the entire (i.e. fitted) column density, we obtain 
a corrected intrinsic flux that would be about 7 times larger (although 
a proper wind absorption model, accounting for the ionization of the 
wind material would be required, see e.g. Waldron et al. 1998). 


(Fig. 4) corresponds to the X-ray emission during the flare 
event, the contribution of the quiescent time intervals to the 
spectrum being essentially negligible. 

3. The nature of the X-ray sources in NGC 6530 

We have plotted the X-ray sources with a single optical 
counterpart in a (B - V,V) colour-magnitude diagram 
(Fig. 5). For this purpose, we adopted < E(B - 17) >= 0.30, 
R v = Av/E(B-V) = 3.1 and DM = 11.25 (van den Ancker 
et al. 1997). Although they found that some highly reddened 
stars show anomalous extinction, van den Ancker et al. (1997) 
proposed that the average extinction law of the intercluster 
material is normal (Rv = 3.1). 

We caution that there could be a number of fore- 
ground (or background) sources that contaminate our sample. 
Nevertheless, for the sake of simplicity, we have adopted the 
distance of NGC 6530 for all the sources. In Table 1 , we include 
the membership probability from the proper motion study of 
van Altena & Jones (1972) whenever this information exists. 
However, Sung et al. (2000) caution that the proper motion 
study might be problematic since low membership probabili- 
ties were assigned to some early-type stars. In their study. Sung 
et al. therefore favour purely photometric criteria. 

Apart from a group of rather bright (V < 12) stars that are 
probably either early-type OB stars, Herbig Ae/Be stars (see 
below) or foreground objects, we find a group of objects ly- 
ing to the right of the ZAMS (by about 0.5 in B - V). Note 
that this result is preserved if we adopt the reddening properties 
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EPIC lightcurve of SCB731 



-1 1 1 1 1 i 1 1 1 L- J 1 I l 1 1 t 1 I i I i i 

0 5 10 15 20 

Time Elapsed [1000s] 


Fig. 3. EPIC light curves of the X-ray source XMMU 
J180435.5— 242729 associated with SCB73I. We have also ex- 
tracted the light curve of a nearby background region (not shown 
here). The latter does not display any variability. 

suggested by Sung et al. (2000) instead. Three of these objects 
display Ha emission and are classified as PMS stars (Sung et 
al. 2000). It seems therefore very likely that this family of X- 
ray sources corresponds to intermediate- and low-mass PMS 
objects. 

3. 1. A population of pre-main sequence stars? 

Low-mass PMS stars are classified according to their infrared 
properties into classes 0, I, II and III (see e.g. Feigelson & 
Montmerle 1999). Class 0 and I objects are deeply embedded 
in their nascent molecular cloud and due to their heavy cir- 
cumstellar extinction, the X-ray emission from these objects is 
strongly attenuated. These protostars have X-ray luminosities 
that range from 10 28 to lO 30 ergs -1 with occassional flares of 
order 10 31 ergs -1 (Carkner et al. 1998) and are best detected 
in the hard X-rays (e.g. Feigelson & Montmerle 1999). 
Class II pre-main sequence stars are surrounded by a thick 
disk. The most distinctive feature of these so-called classical 
TTauri stars is their Ha emission. Sung et al. (2000) used this 
criterion to identify 58 PMS stars in NGC 6530 with strong 
Ha emission and 17 PMS candidates with weak Ha emission. 
As the PMS stars evolve, it is thought that their disks dissipate 
and the PMS stars become Class III objects (or weak-line 
TTauri stars) that have weaker or no Ha emission. At least 


SCB731: EPIC-M0S2 



Energy (keV) 


Fig. 4. EPIC-MOS2 spectrum of the X-ray source XMMU 
J 180435.5-242729 associated with SCB731. The solid line yields 
the best fitting absorbed mekal model. 


in the Taurus- Auriga-Perseus complex, the weak-line TTauri 
stars are found to be X-ray brighter than the classical T Tauri 
stars (Stelzer & Neuhauser 2001). 

The X-ray emission of TTauri stars is probably due to 
solarlike magnetic activity or a magnetic interaction with the 
surrounding protostellar disk. However, the X-ray luminosities 
are up to 10 3 times larger than that of the Sun (Neuhauser 
1997). Low-mass PMS stars in the TTauri stage have multi- 
temperature thermal X-ray spectra with a soft component 
in the range kT t G [0.2, 0.45] keV and a harder component 
with kT 2 € [1.3, 2.6] keV. These objects are usually variable 
X-ray emitters that exhibit flare events with a fast rising curve 
followed by a slower decay (see e.g. Feigelson & Montmerle 
1999 for a review on their high-energy properties). 

In the lower panel of Fig. 5, we have plotted the unab- 
sorbed luminosities (in the 0.5 - 5.0 keV band) inferred from 
the spectral fits in Table 4 and assuming a distance modulus 
of 11.25. For those sources for which we were not able to 
perform a spectral fit, we have converted the observed count 
rates into unabsorbed fluxes assuming a thermal plasma 
model with kT = 2keV, an interstellar column density of 
N h = 5.8 10 21 E(B -V) = 0.17 10 22 cm -2 (Bohlin et al. 
1978) and DM = 11.25. Note that kT = 2keV is in better 
agreement with the mean temperature inferred from our 
spectral fits than the ‘usual’ kT = 1 keV. Adopting kT = 1 keV 
instead of 2 keV would reduce the luminosities by a factor 0 56 
(-0.25dex). 

The unabsorbed X-ray luminosities of the brightest sources are 
found to cluster around a few times 10 31 ergs -1 , as expected 
for quiescent PMS stars (Neuhauser 1997). The strong flare 
that we observe in the light curve of SCB 731 provides further 
evidence that at least some of these sources are related to PMS 
stars. 
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Fig. 5. Top: colour-magnitude diagram of the X-ray sources in the very 
young open cluster NGC 6530. The reddening vector with Rv — 3.1 
is indicated and the solid line shows the ZAMS relation taken from 
Schmidt-Kaler (1982) with a distance modulus DM = 1 1.25 and red- 
dened with E{B -V) = 0.30 (see text). Black dots stand for stars 
from Sung et al. (2000) without Ha emission, asterisks stand for pre- 
main sequence stars or PMS candidates from Sung et al. and open di- 
amonds indicate two stars classified as Be stars by van den Ancker 
et al. (1997). Bottom: unabsorbed X-ray luminosity (in the energy 
range 0.2 - 5.0 keV) as a function of the observed V magnitude. 
The MOS count rates were converted into luminosities assuming a 
2 keV thermal plasma model, a uniform interstellar column density 
of Nh = 0.17 10 22 cm -2 and a distance of 1.8 kpc, except for the 
sources for which a detailed spectral fit is available. For these sources 
(symbols surrounded by an open circle) we used the flux from the 
actual fit. The luminosity of 9 Sgr is taken from Rauw et al. (2002). 
The dashed lines yield different values of log (Lx/Lboi) (indicated 
by the labels) corresponding to main sequence stars of magnitude V 
in NGC 6530. 


In summary, among the objects with an intrinsic luminos- 
ity exceeding 10 31 ergs -1 and having a single optical coun- 
terpart, there are 49 X-ray sources with an optical counterpart 
fainter than 12th magnitude and we suggest that most of these 
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objects are good candidates for TTauri stars. The fact that only 
a few X-ray selected PMS candidates display an Ha emission 
is in agreement with the suggestion that weak-line T Tauri stars 
are intrinsically X-ray brighter than the classical TTauri stars 
(Stelzer & Neuhauser 2001). 

3.2. Early-type stars 

Five Be stars in NGC 6530 were found to exhibit a strong IR 
excess and were suggested to be Herbig Ae/Be stars by van 
den Ancker et al. (1997): LkHa 108 (SCB240), LkHa 112 
(SCB583), LkHa 115 (SCB879), Walker 29 (SCB427) and 
Walker 303. These stars are believed to be intermediate-mass 
(~ 5.5 - 14 Mq, van den Ancker et al.) young stellar objects 
(YSOs). X-ray emission from Herbig Ae/Be YSOs has been 
reported by Hamaguchi et al. (2001). The objects observed by 
Hamaguchi et al. (2001) were found to display high tempera- 
ture thermal emission (kT ~ 1 - 5 keV) with luminosities in 
the range 10 30 - 10 32 ergs -1 . 

Walker 29 and Walker 303 are not detected in our soft- 
band X-ray images. LkHa 112 is marginally detected (with 
a combined likelihood of 9.9), but another star (SCB581) 
lies also inside the 12" radius around the position of the 
XMM source. LkHa 108 lies inside the ‘error box’ of 
XMMUJ180350.5-242110, together with SCB228 and 
229. The only unambiguous detection of X-ray emission 
from a Herbig Ae/Be star in our sample concerns LkHa 1 15 
(XMMUJ 180450.6-242540), which is detected with an 
unabsorbed X-ray luminosity of 3.9 10 31 ergs -1 . We also 
note that the X-ray spectrum of XMMUJ 180345. 1—242205 
(SCB 182, B2.5 V) yields a surprisingly high temperature (kT 
~ 3.2 keV) for a B-star. This suggests that SCB 182 could also 
be an intermediate-mass YSO, despite the fact that van den 
Ancker et al. (1997) did not classify it as a Herbig Be star. 

High-mass young stellar objects could also emit X-rays. 
Chandra observations of the Orion Nebula Cluster revealed 
a huge number of X-ray emitting PMS stars over a broad 
range of masses (Garmire et al. 2000) including the detection 
of emission from massive stars that are just settling down on 
the ZAMS. Rho et al. (2001) reported on ROSAT and ASCA 
observations of the Trifid Nebula (M20) which is thought to 
be in a ‘pre-Orion’ star forming stage. Two massive YSOs as 
well as one low-mass PMS star were identified as counterparts 
of ROSAT sources. Kohno et al. (2002) obtained a Chandra 
observation of the MonocerosR2 cloud and detected X-ray 
emission from four massive YSOs which were found to have 
a mean temperature of kT ~ 3 keV and to display flare-like 
events, much like their low-mass counterparts. This led Kohno 
et al. (2002) to suggest that the massive YSOs produce X-rays 
through the same magnetic activity as low-mass PMS stars. 
This activity would continue until the ZAMS phase and the 
onset of the stellar wind. 

In NGC 6530 it seems rather unlikely to find high-mass YSOs 
that have not yet reached the ZAMS (at least outside the HG 
region). In fact, the earliest star in our field of view (apart 
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from 9Sgr and H36) is HD 164816 (09.5 IIMV) which has 
probably already evolved off the main sequence. Accordingly, 
van den Ancker et al. (1997) suggested that the formation 
of stars more massive than 8M 0 has already stopped in 
NGC6530. An exception to this rule could be the Hourglass 
Region which we shall discuss in the next section. 

Finally, it is quite remarkable that some of the optically 
brightest (V < 10.0) stars in NGC6530 are not detected in 
our EPIC data. These include the late-type stars HD 164584, 
HD 164948 and SCB444, as well as a number of main- 
sequence B-type stars (HD 164865, HD 164947, SCB 466, 588, 
599, 647, 667 and 708) that have spectral types B1 or slightly 
later according to van den Ancker et al. (1997). The latter result 
is somewhat surprising. In fact, ROSAT all-sky survey results 
reported by Berghofer et al. (1997) suggest that stars of spec- 
tral type B1 can have X-ray luminosities of order a few times 
10 31 ergs” 1 and most of them should therefore be detected in 
our data. We conclude that the early B-type main sequence stars 
in NGC 6530 may have X-ray luminosities on the lower side of 
the Lx/Lboi relation proposed by Berghofer et al. (1997). 

4. The Hourglass Region 

Our EPIC data also unveil X-ray emission from the Hourglass 
Region in M 8. Beside the ionizing star Herschel 36, the HG 
also harbors the ultra-compact Hll region G5.97-l.17 as 
well as a number of infrared sources some of which are 
believed to correspond to embedded massive stars. In fact, 
Woodward et al. (1990) performed near infrared imaging of 
the HG, unveiling a Trapezium-like stellar cluster of very 
young hot stars at the core of the HG. G 5. 97 -1.17 could 
be a young stellar object (probably a late-B star) surrounded 
by a circumstellar disk that suffers photoevaporation by the 
radiation from Herschel 36, similar to the so-called proplyds 
in the Orion Nebula (Stecklum et al. 1998). Herschel 36 is 
probably a massive YSO that is more evolved than the exciting 
stars of ultra-compact Hll regions. In fact, Stecklum et al. 
(1995) reported infrared and HST Ha images that reveal an 
elongated jet-like structure associated with H36 orientated 
roughly perpendicular to the elongated dust structure around 
the star that could be indicative of a circumstellar disk. 

Woodward et al. (1986) proposed a ‘blister’ scenario of the 
star formation in the HG where the massive star H 36 formed 
near the edge of the molecular cloud. This picture is also 
supported by the results of Chakraborty & Anandarao (1997, 
1999). These authors studied the kinematics of the HG and 
found the H II bubble around H 36 to be expanding which sug- 
gests that the ionization front is still destroying the surround- 
ing molecular cloud and pressure equilibrium has not yet been 
reached. In addition, these authors reported high velocity flows 
to the south of H 36 that indicate a Champagne flow. 

The EPIC soft-band image of the HG region (Fig. 6) 
reveals an extended X-ray emission probably due to the 
contributions of the young stellar object H36, the UC Hll 
region and possibly a diffuse emission from a hot bubble 


the Lagoon Nebula and the veiy young open duster NGC 6530 



Right Ascension (J2000) 

Fig. 6. EPIC X-ray contours superimposed on a HST WFPC2 Ha 
image of the HG region in the Lagoon Nebula. The X-ray contours 
reveal an elongated X-ray source at the location of the 07 V star 
Herschel 36 and the UC Hll region G 5.97 -1.17 as well as an emis- 
sion lobe extending over the southern part of the optical HG nebula. 
The EPIC contours were smoothed by convolution with a gaussian of 
a = 2 pixels (1 pixel = 1 arcsec 2 ) and the contour levels correspond 
to 0.3, 0.4, 0.55, 0.81 and 1 EPIC (MOS1 + MOS2 + pn) count per 
arcsec 2 in the 0.5 - 1.2 keV band. The X-ray source below the HG is 
associated with SCB 142 and 146 (see Table 2). 


created by the wind of Herschel 36. In particular, we find no 
correlation between the soft X-ray emission and the embedded 
IR sources discussed by Woodward et al. (1990). While the 
X-ray emission peaks around the position of H 36, it presents 
also an enhancement towards the south-east of H36 (roughly 
coincident with the southern lobe of the optical HG nebula). 
Woodward et al. (1986) found that the electron density 
(4000 cm” 3 ) on the eastern side of the ionized cavity is about 
ten times larger than the electron density on the western side 
(see their Fig. 15). Therefore, we find that the X-ray emission 
is enhanced over the higher density part of the HG region. 

We have extracted EPIC-MOS spectra of the Hourglass in 
the energy range 0.2 - 1.2 keV and fitted them simultaneously 
with the ROSAT-PSPC spectrum extracted from observation 
rp900374n00. The source counts were extracted over an 
elliptical area designed to avoid contamination from neigh- 
bouring sources. The background-corrected spectra can be 
fitted (xl = 1.25, 55 d.o.f.) with an absorbed mekal model 
with N h = (1.11±;}*) 10 22 cm” 2 , kT= 0.631 .05 keV. The ob- 
served and unabsorbed fluxes in the energy band 0.2 - 2.0 keV 
are respectively 1.1 10” 13 and 1.7 10” 12 ergcm~ 2 s” 1 . The 
latter value corresponds to an integrated intrinsic luminosity of 
6.6 10 32 erg s" 1 . Assuming that the X-ray emission from H 36 
roughly follows the Lx/Lboi relation proposed by Berghofer et 
al. (1997), we expect H36 to contribute 1.5 lO^ergs” 1 , i.e. 
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about one fourth of the total X-ray luminosity of the HG region. 

Although we cannot completely rule out that the shape of 
the contours in Fig. 6 might result from the superposition of 
the PSFs of several faint point sources, it seems very tempting 
to associate the soft X-ray emission of the HG with a diffuse 
emission from a hot plasma around H 36. Diffuse X-ray 
emission has been observed in a few Hu regions only. For 
instance, Wang (1999) reported diffuse X-ray emission arising 
in blister shaped regions of the 30 Dor complex outlined by 
loops of Hu gas. Such a description could also suit the X-ray 
emission from the Hourglass Region. On the other hand, 
recent Chandra observations of the massive galactic starburst 
NGC 3603 revealed an extended diffuse X-ray emission in 
the immediate surroundings of the cluster core (Moffat et 
al. 2002). Moffat et al. (2002) suggest that the diffuse X-ray 
emission may arise from multiple merging/colliding hot stellar 
winds, though the exact contribution of unresolved PMS stars 
is not fully constrained. Since the HG region is a much less 
extreme star formation region, a superwind scenario (as in 
NGC 3603) that would require a large concentration of hot 
massive stars, is less likely to account for the X-ray emission 
seen in the HG. 

As pointed out above, the density distribution in the 
vicinity of the Hourglass is believed to be highly non-uniform 
(Woodward et al. 1986) and a sophisticated numerical model 
such as the one developed by Comeron (1997) is therefore 
needed for a detailed comparison. In fact, simulated X-ray 
maps of a blister model presented by Comeron (1997) show 
that the X-ray emission consists of two components: an 
extended one from the bubble expanding into the intercloud 
(low-density) medium and a compact one in the vicinity of the 
star (i.e. in the high density part still contained in the molecular 
cloud). The simulations of Comeron indicate that the intrinsic 
soft X-ray emission from the compact component should have 
an emission peak approximately 40 times more intense than 
the extended emission component. The observed enhance- 
ment of the extended emission over the high density part of 
the HG is at least in qualitative agreement with this description. 

Could the stellar wind of Herschel36 (07 V) account for 
the X-ray luminosity of a wind-blown bubble? The theory 
of these bubbles has been considered by Castor et al. (1975) 
and Weaver et al. (1977), with further additions by Mac Low 
& McCray (1988) and Chu et al. (1995) amongst others. 
The stellar wind sweeps up a dense shell of material, but 
within the shell there will be a region of hot shocked gas 
(see e.g. Strickland & Stevens 1999 for a schematic view 
of the structure of a wind-blown bubble). Within the frame- 
work of this theory, the parameters of the bubble depend 
upon the wind luminosity of the star, the density of the ambi- 
ent medium and the age of the bubble (see e.g. Chu et al. 1995). 

To derive an order of magnitude estimate, we adopt 
T eff = 39000 K and log(L bo i/L 0 ) = 5.31 for H36. 
From these parameters, we expect a mass-loss rate 
M = 4 1O~ 7 M 0 yr -1 and a terminal wind velocity 


Voo = 2300 km s -1 (Howarth & Prinja 1989, Prinja et 
al. 1990). This yields a kinetic luminosity of the star 
(£wind = l/2Mi^o) of 6.7 10 35 ergs -1 . For simplicity, we 
assume that the ambient density is constant at the higher value 
derived by Woodward et al. (1986) with n e = 4000 cm~ 3 . 
Adopting an expansion velocity of the HG of ~ lOkms -1 
(Chakraborty & Anandarao 1997, 1999) and an angular radius 
of the cavity in the molecular cloud of about 25" (correspond- 
ing to 0.2 pc at a distance of 1.78 kpc), we infer a dynamical 
age of about 20000 yrs. 

From these values, we derive the theoretical properties of 
the bubble surrounding H 36. The central temperature of the 
bubble should be ~ 5 10 6 K (i.e. kT ~ 0.44 keV), the predicted 
radius is 0.4 pc and the expected X-ray luminosity amounts to 
Lx ~ 10 35 ergs _1 . 

We note that the predicted radius is larger than the observed 
radius. However, the discrepancy is less pronounced than for 
some other cases in the literature (e.g. Naze et al. 2001 and 
references therein). The most outstanding discrepancy between 
the model and our observations of the HG region concerns the 
expected X-ray luminosity which is more than a factor 175 
larger than the observed one. Part of this discrepancy may be 
due to our estimate of the kinetic luminosity of the wind of 
H 36. For instance, Stecklum et al. (1998) noted that the WE 
spectra of H 36 are similar to those of 6 1 Ori C (HD 37022) 
and suggested hence that H 36 is close to the ZAMS and 
has a less developed wind than a typical 07 main sequence 
star. Assuming a mass-loss rate four times lower and a wind 
velocity half the values used above, the discrepancy would be 
reduced to a factor rKj 12. Note however that Leitherer (1988) 
derived M = 7.6 10 -7 M© yr _1 and Voo = 1650 km s -1 for 
6 1 OriC which yields about the same value of the kinetic 
luminosity than derived hereabove. Another parameter that 
might account for some of the discrepancy is the density of the 
ambient medium. We have adopted a rather large value which 
should be appropriate for the molecular cloud. However, in 
the blister scenario, part of the bubble expands into the much 
lower density intercloud medium and the X-ray luminosity of 
the bubble around H 36 may be set by the density of the latter. 

Although the theoretical values are way above the observed 
X-ray luminosity, we conclude that the kinetic luminosity of 
the stellar wind of Herschel 36 is most probably sufficient to 
account for the existence of a wind-blown bubble. 

5. Summary and conclusions 

Our XMM-Newton observation of the Lagoon Nebula (M 8) 
reveals a cluster of rather faint point-like X-ray sources 
in addition to the bright 04 star 9Sgr. We have shown 
that most of these sources have optical counterparts inside 
the very young open cluster NGC 6530. We suggest that 
an important fraction of these sources are associated with 
low-mass pre-main sequence stars (most probably weak-line 
TTauri stars) which were not previously identified as PMS 
objects because of the lack of strong Ha emission. Sung 
et al. (2000) estimated that the detection efficiency from a 
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single X-ray observation amounts to ~ 46%, quite similar 
to that of a single Ha observation. Our results suggest that 
the true number of low-mass PMS objects exceeds by more 
than a factor two the number of objects identified by Sung et al. 

The X-ray selected PMS candidates in NGC 6530 are 
found to have rather hard spectra with kT of a few keV and 
their X-ray luminosities (of order of a few times 10 31 erg s' 1 ) 
are apparently in the higher part of the luminosities of TTauri 
stars in the literature. The large number of X-ray sources 
associated with PMS stars in NGC 6530 may shed some light 
on the origin of apparently diffuse emission in more distant 
starburst clusters (e.g. NGC 3603) where source confusion 
becomes a serious issue. 

As to intermediate and high-mass young stellar objects, 
we found that at least one of the known Herbig Be stars in 
NGC 6530 exhibits a relatively strong X-ray emission, while 
most of the early-type main sequence stars of spectral type B1 
and later are not detected. 

We also report X-ray emission from the massive proto- 
star Herschel 36 as well as probably diffuse emission from the 
Hourglass Region that might reveal a bubble of hot gas pro- 
duced by the interaction of the stellar wind of Herschel 36 with 
the denser part of the molecular cloud. Higher spatial resolu- 
tion observations of the HG with Chandra are needed to con- 
firm the diffuse nature of this emission and to further constrain 
its morphology. 
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